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Abstract 

Exclusive semileptonic B decays into excited charmed mesons are studied 
with QCD sum rules in the leading order of heavy quark effective theory. 
Two universal Isgur-Wise functions r and £ for semileptonic B decays into 
four lowest lying excited D mesons (D±, D\, D' , and D[) are determined. 
The decay rates and branching ratios for these processes are calculated. 
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I. INTRODUCTION 



The heavy quark effective theory (HQET) |],|2| is a useful tool to describe the spec- 
troscopy and matrix elements of mesons containing a heavy quark. In the infinite mass 
limit, the spin and parity of the heavy quark and that of the light degrees of freedom 
are separately conserved. Coupling the spin of light degrees of freedom ji with the spin 
of heavy quark sq = 1/2 yields a doublet of meson states with total spin j = ji ± 1/2. 
The ground state mesons with Q q flavor quantum numbers contain light degrees of free- 
dom with spin-parity j[ = | , giving a doublet (0~,1 - ). For Q = c these mesons are 
the doublet (D,D*), while Q = b gives the doublet (B,B*). The excited charmed mesons 
with j[ = l + /2 and 3 + /2 are two spin symmetry doublets (0 + ,l + ) and (l + ,2 + ). These 
mesons are doublets (D' ,D[) and (D^D^)- The D\ and D\ mesons have been observed 
with rather small widths, while D' and D[ have not been observed for the reason of being 
too broad. The properties of these lowest lying excited mesons have attracted attention 
in recent years. The mass and decay widths have been studied with potential model ]3].Q . 
relativistic Bethe-Salpeter equation f5j and QCD sum rules [| f| . 

Progress has been achieved recently in the studies of semileptonic B decays into excited 
charmed mesons {D' Q , D[, D\ and D^)- Semileptonic B decay into an excited heavy meson 
has been observed in experiments P Hl^l - With some assumptions, CLEO and ALEPH 
collaborations have reported respectively the branching ratios BiB — > Dieu e ) = (0.56 ± 
0.13 ± 0.08 ± 0.04)% and B{B -»■ D x ev e ) = (0.74 ± 0.16)%, as well as the limits B{B -»■ 



B\ev e ) < 0.8%(90% C.L) and B(B -»■ D*eP e ) < 0.2% |K|[12|. Theoretically, HQET 
provides a systematic method for investigating such processes. The semileptonic B decay 
rate to an excited charmed meson is determined by the corresponding matrix elements of 
the weak axial-vector and vector currents. Heavy quark symmetries can be used to reduce 
the form factors parameterizing the matrix elements. In the infinite quark mass limit 
these matrix elements are described respectively by one universal Isgur-Wise function and 
vanish at zero recoil |13|,I^|. Extensive investigation in |14| shows that the leading 1/mg 
correction at zero recoil can be calculated in a model independent way in terms of the 
masses of charmed meson states. 

The universal function embodies details of low energy strong interactions and cannot be 
calculated from first principles. It must be calculated in some nonperturbative approaches. 
For that purpose, there are many viable approaches, including different quark models 
15H21|, relativistic Bethe-Salpeter equation |22[ and QCD sum rules |23Vp6|. In the present 



work we study the B semileptonic decays to excited charmed meson states (D' ,D[) and 
(DijD^) with QCD sum rule in the leading order of HQET. In particular we compute the 
relevant universal Isgur-Wise functions that describe such decays in the mq — > oo limit. 
The remainder of this paper is organized as follows. In Section [TI| we present the 
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formulae of weak current matrix elements and of decay rates. In Section [II we begin with 
a brief review on the interpolating currents for excited heavy mesons and on relevant sum 
rules for two-point correlators, then the sum rules for Isgur-Wise functions r and ( are 
derived. Section [TV] is devoted to numerical results and implications. Finally in Section |TV| 
we draw our conclusions. 



II. ANALYTIC FORMULAE FOR SEMILEPTONIC DECAY AMPLITUDES 

B -> (D u D* 2 )lv AND B -> (D' ,D[)£u 

The theoretical description of semileptonic decays involves the matrix elements of vector 
and axial vector currents (V M = c^h and = 07^756) between B mesons and excited 
D mesons. For the processes B — > Diiu and B — ► D\lv, these matrix elements can be 
parameterized as 

(D 1 (v f , e)| V* \B(v)) = yJm D , m B [f Vl ^ + (fay* + fvj») e* ■ v] , (la) 
(Dxiv', e)| A" |B(v)> = i m B / A e^e>; , (lb) 

e)| A" |B(«)> = 7m D . m B e*^ a + (k A j>* + v a v p ] , (lc) 

(£>£(t/,e)|V|B(t;)) =«7m^%fcyrt>V;. (Id) 

The form factors /j and fcj are functions of y — v ■ v'. In the limit tuq — ► 00 they can be 



expressed by a universal dimensionless function r(y) |T3|,|T4 
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Mv) = 77= (1 - V 2 ) r{y) , = --^r(y), 

A*(v) = 4 (y - 2 ) , = (1 + v) r(») , 

^(y) = -(l + 2/)r(y) , fc A2 (y)=0, 

fc A 3 (l/) = , k v (y) = -r(w). (2) 

Since the polarisation vector of the spin-one state D\ and the polarization tensor of the 
spin-two state D\ satisfy e* • v' = and e* m v' a = respectively, only the form factor f Vl 
contributes at zero recoil (v = v'). 

The form factors that parameterize the B — > D'^lv and B — > D' x iv matrix elements of 
the weak currents are defined by 

(D> o (v')\V»\B(v)) = 0, (3a) 

(D' (v')\ A» \B(v)) = ^/m D[ m B [g + {v» + i;'") + g_ («" - «'")] , (3b) 

(D[(v', e)\ V \B(v)) = ^m D[ m B [g Vl e*^ + {g v ^ + g Vz v'») (e* ■ v)] , (3c) 

(D[(v', e)\A"\B(v)) = y/m^rriB % g A e> ^ , (3d) 
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where g$ are functions of y. At zero recoil the matrix elements are determined by g+(l) 
and pvi(l). In the infinite mass limit form factors gi can be written in terms of a single 
function ((y) |TjyTJ], 



g+(y) = o, g-(y) = ((y), 
9A(y) = ((y), 9vAy) = (y- l)C(y) 
»«a(y) = 0, g V3 (y) = -t(w) . 



(4) 



Note that the notations of Ref. |14[] have been used here, r is v^3 times the function r 3 / 2 
of Ref. ||13|| , while ( is two times the function of Ti/ 2 . All of these matrix elements of the 
weak currents vanish at zero recoil y = 1, since the B meson and the (Di, or (D' Q , D[) 
mesons are in different heavy quark spin symmetry multiplets, and the current at zero 
recoil is related to the conserved charges of heavy quark spin-flavor symmetry. 

The differential decay rates are given by (taking the mass of the final lepton to zero) 



dr 



dy 



G2 ^ Vcbl2m ^rf(y + l)(y 2 -lf 2 



72tt 3 



(y- + n) 2 + y(l - 2r lV + r 2 ] 



dT D * G 2 F \V cb \ 2 m% 3 



dy 



72tt 3 



4 (y + \){y 2 - if' 2 (y + !)(!- r 2 ) 2 + 3y(l - 2r 2 y + r 2 2 ) 



r(y)\\ 
\r(y)\ 2 



dT D , _ G 2 F \V cb \ 2 m% 

dy 487r 3 
dlV G 2 F \V cb \ 2 m% 



r' 3 (y' 



\3/2 



^2 



IC(y)l 5 



r[ 3 (y-l)Jy 2 -l (y - 1)(1 + r[) 2 + %(1 - 2r[y + r 2 ) \C{y)\ 



'2\ 



dy 487r 3 

where 7*1 = m^/ms, r 2 = mD*/mB, r' Q = m D JmB and r[ = m^/ms- 



(5) 
(6) 
(7) 
(8) 



III. SUM RULES FOR ISGUR-WISE FUNCTIONS r AND ( 

A. Interpolating currents for heavy mesons of arbitrary spin and parity and 

two-point correlation function 

A basic element in the application of QCD sum rules to excited heavy mesons is to 
choose a set of appropriate interpolating currents in terms of quark fields each of which 
creates (annihilate) a definite excited state of the heavy mesons. The proper interpolating 
current J^p j/ for the state with the quantum number j, P, ji in HQET was given in 0. 
These currents have nice properties. They were proved to satisfy the following conditions 

(0 ./;--(())/./>'../;) =ifp Sl 5 ji ,8pp,5 j j7f*™°' , (9) 
i ( l T (^r^)^,f ''(°)) 1°) = S^pp^i-iy S ■ ■ ■ g^ 

x / dt6(x-vt)Up tk (x) (10) 
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in the tuq — ► oo limit. Where rj ai "' a j is the polarization tensor for the spin j state, v is the 
velocity of the heavy quark, g^ 13 = g af3 — v a v 13 is the transverse metric tensor, S denotes 
symmetrizing the indices and subtracting the trace terms separately in the sets (ai • • • atj) 
and (fli • • • j3j), fpj e and Hpj e are a constant and a function of x respectively which depend 
only on P and j^. Because of equations @ and (|IUD, the sum rules in HQET for decay 
amplitudes derived from a correlator containing such currents receive contribution only 
from one of the two states with the same spin-parity (j, P) in the tuq — > oo. Starting from 
the calculations in the leading order, the decay amplitudes for finite tuq can be calculated 
unambiguously order by order in the I/uiq expansion in HQET. 

In the following we focus our attention on the semileptonic decays, B — > Di, D% and 
B — > D' Q , D[. The relevant doublets to be considered are ground states and lowest lying 
positive parity states, namely, doublets (0~,1~), (0 + ,l + ) and (l + ,2 + ). The currents for 
creating 0~ and 1~ are usual pseudoscalar and vector currents 

Jl a 1 = \ - Ki^q , J] a i = \l~ Kn?Q ■ I I 




As pointed out in 0, there are two possible choices for currents creating + and 1 + of the 
doublet (0+,l+), either 

4,+,2 = h v q , (12) 

J\X, = ^ Kl'^q , (13) 

or 

$,+,2 = ^=hv{-i)PtQ, (14) 

j'i] + ,2 = A= h vl ^{-i)lp t q , (15) 

where D is the covariant derivative. Similarly, there are two possible choices for the currents 
creating 1 + and 2 + of the doublet (l + ,2 + ). One is 

= \/f Kl\~i) {V? - ^Pt) Q , (16) 




4+;r = \ll n v [ -T- (V^ 2 + ^ v * x - \ 9t ia2 ft) q ■ (1 7 ) 



Another choice is obtained by adding a factor — %p t to flT6|) and (p~T|) . Note that, without 
the last term in the bracket in (p!6|) the current would couple also to the 1 + state in the 
doublet (0 + , 1 + ) even in the limit of infinite tuq. 
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Usually, the currents with the least number of derivatives are used in the QCD sum 
rule approach. However, there is certain motivation for using the currents (|HP, (O) for 
the doublet (0 + , 1 + ) for the reason that in the non-relativistic quark model, the doublets 
(0 + , 1 + ) are orbital p-wave states which correspond to one derivative in the space wave 
functions. Therefore, we shall consider both the currents (|12"D, ( pT3[ ) and (0), flTop for the 
doublet (0+ 1+). 

For the ground state heavy mesons, the sum rule for the correlator of two heavy-light 
currents is well-known. It is: [p|,[25[| 



f. i e - 2A - /T = 4- 2 r - 2 ^ /T d- - \m f i - m 



16?r 2 Jo 2^' \ AT 2 J v ' 

For the doublet (0 + , 1 + ), when the currents J'o,+,2, J'1,+,2 in (0), (|15|) are used the sum 
rule (same for the two states) after the Borel transformation is found to be @ 

The corresponding formula when the current Jq,+,2 and Ji,+,2 in (|12|) and flIE|) are used 
instead of J'0,+,2 and J'i,+,2 is the following 

/^•T^ = JL f Je-^d. + i <»> - mS (ft) . (20) 

When the currents (|IED and ( fPT| ) are used the sum rule for the 2 + ) doublet is found 
to be 

fl 3e" 2A -i /T = -±* r ^e-^dco -^ml (qq) - 4^G 2 )T . (21) 

J+,- 2 2 6 7T 2 Jo 12 ° Wy/ 2 5 V ' K ' 

Here Wq (gg) = (qga^G^q). Above results will be used in the following sections. 



B. sum rules for r and £ 



For the ampltudes of the semileptonic decays to excited states in the infinite mass limit, 
the only unknown quantities in the expressions (5)-(8) are the universal functions r(y) and 
C(y)- Let us first consider the Isgur-Wise function r. In order to calculate this form factor 
using QCD sum rules, one studies the analytic properties of the three-point correlators 

i 2 j fxfa (0|T (j v lM {x) J^\0) J{_i{z)) |0> = E(cu,u;',y) C% A , (22a) 

i 2 J d'xd'ze^-^ <0|T (j£s(x) ^ ( r°(0) |0> = E{u,u/,y) £%$ , (22b) 

where Jy^ v ' v ^ = h{y')^ h(v), J^} v,v ■* = h(v') r y fl 'y 5 h(v). The variables k, k' denote residual 
"off-shell" momenta which are related to the momenta P of the heavy quark in the initial 
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state and P' in the final state by k = P — rriQV, kl = P' — txiqiv' respectively. For 
heavy quarks in bound states they are typically of order Aqcd and remain finite in the 
heavy quark limit. Cy.A are Lorentz structures associated with the vector and axial vector 
currents(see Appendix). 

The coefficient H(u;,a/, y) in (|22| ) is an analytic function in the "off-shell energies" 
uj = 2v ■ k and uj' = 2v' ■ k' with discontinuities for positive values of these variables. It 
furthermore depends on the velocity transfer y — v ■ v', which is fixed at its physical region 



for the process under consideration. By saturating (|22| ) with physical intermediate states 
in HQET, one finds the hadronic representation of the correlator as following 
_ f-if + s.r(y) 

Z hadro (uj, uj', y) = —j ' 2 . : r-r + higher resonances , (23) 

(2A_ ) i — uj — «e)(2A +j | — ur — it) 

where fpj e are constants defined in (|9]), Apj e = m P j e — uiq. As the result of equation (|j), 
only one state with j p = 1 + or j p = 2 + contributes to (|23|), the other resonance with the 
same quantum number j p and different ji does not contribute. This would not be true for 
j p = 1 + if the last term in (0) is absent. 

Using the expression for the heavy quark propagator / dt5(x — vt)(l + f)/2 in HQET 
the correlators in (^) have the following form 

1 + j> l + fl 



J dhdt 2 e ^ + ^ 2 Tr -^z 7 ^^l(o|r' g (^2)g(-^i)r|o> 



where T and T' are functions of Dirac matrices and covariant derivatives appearing in the 
definations of the currents. For u and uj' in the deep Euclidian region the integral is dom- 
inated by small values of t\ and t 2 . As usually done in QCD sum rule approach, uj and 
uj' are analytically continued to the deep Euclidian region. Therefore OPE in the short 
distance can be applied to the above expression for any value of y not much larger than 1. 
This is because, for small tx, t 2 and y of the order 0(1), all components of vt\ and v't 2 can 
be made small by choosing the rest frame of B. In our case the maxima value of y is less 
than 1.5. Therefore S(uj,uj',y) in ( [22] ) can be approximated by a perturbative calculation 
supplemented by nonperturbative power corrections proportional to the vacuum conden- 
sates which are treated as phenomenological parameters. The perturbative contribution 
can be represented by a double dispersion integral in uj and uj' plus possible subtraction 
terms. We find that the theoretical expression for the correlator has the form: 

Etheo(^,^',y) — / dvdv -. — 7 — --(-subtractions 

J (y — uj — it) [y' — uj' — it) 

+ E cond (u,uj',y) . (24) 

Following the standard QCD sum rule procedure the calculations of E(u>,uj',y) are 
straitforward. Confining us to the leading order of perturbation and the operators with 
dimension D < 5 in OPE, the relevant Feynman diagrams are shown in Fig 1. 
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The perturbative part of the spectral density is 

3 1 / 

p pert (u,u',y) = — — ( - 3^ 2 + (^' 2 + 2uu'){2y - 1)) 

x G(w) G(cD') G(2yuu' - £ 2 - cD' 2 ) . (25) 

The QCD sum rule is obtained by equating the phenomenological and theoretical ex- 
pressions for S. In doing this the quark-hadron duality needs to be assumed to model the 
contributions of higher resonance part of Eq. fl23|) . Generally speaking, the duality is to 
simulate the resonance contribution by the perturbative part above some threshold ener- 
gys. In the QCD sum rule analysis for B semileptonic decays into ground state D mesons, 
it is argued by Neubert and Shifman in j|2|,pi],^5|l that the perturbative and the hadronic 
spectral densities can not be locally dual to each other, the necessary way to restore dual- 
ity is to integrate the spectral densities over the "off-diagonal" variable £>_ = (u) — &')/2, 
keeping the "diagonal" variable uj + = (u + <D')/2 fixed. It is in uj + that the quark-hadron 
duality is assumed for the integrated spectral densities. We shall use the same prescription 
in the case of B semileptonic decays into excited state D mesons. 

The G functions in (p5[) imply that in terms of uj + and uj- the double discontinuities of 
the corrrelator are confined to the region —\Jy l — 1/ (1+y) < < Vv 2 ~~ V (1 + 1/) 
and Cj + > 0. According to our prescription an isosceles triangle with the base oj + = uj c is 
retained in the integation domain of the perturbative term in the sum rule. 

In view of the asymmetry of the problem at hand with respect to the initial and final 
states one may attempt to use an asymmetric triangle in the perturbative integral. How- 
ever, in that case the factor (y 2 — l) 3 / 2 in the denominator of (p5|) is not canceled after the 
integration so that the Isgur-Wise function or it's derivative will be divergent at y = 1. 
Similar situation occurs for the sum rule of the Isgur-Wise function for tansition between 
ground states if a different domain is taken in the perturbative integal |25|j . 

In order to supress the contributions of higher resonance states a double Borel trans- 
formation in u and uj' is performed to both sides of the sum rule, which introduces two 
Borel parameters T\ and T 2 . For simplicity we shall take the two Borel parameters equal: 
T\ = T 2 = 2T. In the following section we shall estimate the changes in the sum rules in 
the case of T\ ^ T 2 . 

After adding the non-perturbative part and making the double Borel transformation 
one obtains the sum rule for r as follows 

- -±- (?LGG)y±*. (26) 
3x2 5 ^tt (y + 1) 2 

We now turn to the study of ((y). To obtain QCD sum rules for this Isgur-Wise function 
one starts from three-point correlators 
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i 2 



J d^z^) (Q\T (j , + i(x) jf' v '\0) |0> = H( W ,a/,y) ^ , (27a) 

| d*xd*ze« k '- x - k ^ (0|T (jr, + ,i(^) ^ y) (0) |0> = S^V.y) , (27b) 



the inserting of vector current in ( |27a| ) vanishes. 

As mentioned above there are two possible choices for the currents creating + and 
1 + of the doublet (0 + , 1 + ). Let us first consider the sum rule by using the interpolating 
currents (|TJ|) and flip] ) with derivative. The perturbative part of the spectra density is 
found to be 



p' pert (Co, Co', y) = ^ , 2 _ 1)3/ 2 (ft - ©(«0 0(2^' " & - w ') > 

(28) 

With the same procedure the resulting sum rule for ( takes the form 

cw 4i *- (a -^ )/t - ^ r — ^ - t^ i -¥ 

i { ^ GG> Zv + i (29) 



3 x 2 6 N vr ' y + 1 

When the currents QT2] ) and (|T3p without derivative are used, there is no diagram shown 
in Fig (c). The evaluation of the perturbative graph gives the spectral density 

3 1 j 

Ppert(a), u>', y) = - /2 (Co - u) ®{u) &{&') Q(2yCjCu' - Cj 2 - Co' 2 ) . . (30) 

After rewriting the spectral function in terms of u±, performing the double Borel trans- 
formation and the integral over lj_ in the confined region, we find that the perturbative 
contribution vanishes for 7\ = T 2 . Therefore, the sum rule in this case is not a good one. 
We shall not use it in the numerical analysis in the following. 

We end this subsection by noting that the QCD 0(a s ) corrections have not been in- 
cluded in the sum rule calculations. However, the Isgur-Wise function obtained from the 
QCD sum rule actually is a ratio of the three-point correlator to the two-point correlator 
results. While both of these correlators subject to large perturbative QCD corrections, 
it is expected that their ratio is not much by these corrections significantly because of 
cancelation. This has been proved to be true in the analysis for B semileptonic decay to 
ground state heavy mesons . 



IV. NUMERICAL RESULTS AND IMPLICATIONS 

We now turn to the numerical evaluation of these sum rules and implications. For the 
QCD parameters entering the theoretical expressions, we take the standard values 
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(qq) = -(0.23) 3 GeV 3 , 
(a s GG) = (0.04) GeV 4 , 

ml = (0.8) GeV 2 , (31) 



In the numerical calculations we use the physical masses for B, D\ and D\ [p7 |. Whereas 
for D' and D[ we take the theoretical mass values of the doublet in the leading order 
obtained in 0. Therefore we use 

m B = 5.279, m Dl = 2.422, m D * = 2.459, m D * = m D ' x = 2.4, (32) 

as well as = 0.04. 

In order to obtain information for r{y) and ((y) from the sum rules which is independent 
of specific input values of /'s and A's, we adopt the strategy to evaluate the sum rule by 
eliminating the explicit dependence on the parameters /'s and A's by using the sum rules 
for the correlator of two heavy-light currents. Dividing the three point sum rules in (|26"D 
and Q29] ) by the square roots of relevant two point sum rules in QT%D, flT9D and (PH), we 
obtain expressions for the r and ( as functions of the Borel parameter T and the continuum 
thresholds. This procedure may reduce the systematic uncertainties in the calculation. 

Let us evaluate numerically the sum rules for r(y) and C{u)- Imposing usual criterium 
for the upper and lower bounds of the Borel parameter, we found they have common 
sum rule "window": 0.7 < T < 1.1, which overlaps with those of two-point sum rules 
p5]]7||. Notice that in the case of transiton between ground states the normlization of 
the Isgur-Wise function at zero recoil leads to the requirement that the Borel parameter 
T\ = T 2 = 2T. That is, the Borel parameter in the sum rule for three-point correlators is 
twice the Borel parameter in the sum rule for the two-point correlators. In Fig. 2, we show 
the range of predictions for r(y) and ((y) obtained by varying the continuum thresholds 
lu c , u c0 , u' cl and uj c2 , and fixing the Borel parameter at T = 0.9. In the evaluation we have 
taken 2.0 < lu c < 2.5, 1.9 < lu c0 < 2.4, 2.3 < u' cl < 2.8 and 2.5 < lu c2 < 3.0. Where u c0 , 
uj' c1 and u C 2 are relevant two-point sum rule continuum thresholds defined in ([18]), (|i~9D and 
(pl|), their ranges are determined by two-point sum rule analysises p5,7|. The results are 
stable against reasonable variations of these parameters. 

The resulting curves for r(y) and £(y) may be well parameterized by the linear approx- 
imations 

r(y) = r(l)(l-p 2 T (y-l)), r(l) = 0.74 ± 0.15 , p\ = 0.90 ± 0.05 , (33) 
C(y) = C(l)(l-p 2 (y-l)), C(l)= 0.26 ±0.08, p\ = 0.50 ± 0.05 . (34) 

The errors reflect the uncertainty due to cj's and T. 

We have estimated the variations of Isgur-Wise functions at zero recoil in the case of 
T\ 7^ T2. After Borel transformation the weight function in the dispersion integrals is 
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exp{-{l/T x + l/T 2 )u + - (1/Ti - l/T 2 )cD_}. We found that for T 2 = 37\/2 the value of 
r(l) increases about 10%, while the value of £(1) increases about 30%. The variations are 
slower for 7\ larger than T 2 . This gives us the uncertainties of the results within reasonable 
region of Ti/T 2 . 

The numerical values of r and ( at zero recoil are compared with other approaches in 
Table |. For r we find a broad agreement with some of the constituent quark model results, 
whereas for ( we only agree with | ]21{| . 

Using the forms of linear approximations for r(y) and ((y) we can compute the total 
semileptonic rates and branching rations. The maximal values of y in the four cases are 
Vnax = (1 + rj)/2n « 1.32, = (1 + r 2 2 )/2r 2 « 1.31 and y™L = = (1 + r^ 2 )^ « 
1.33, respectively. The results turn out to be 

r Dl = 1.36 x HT 15 GeV , 

T D2 = 2.11 x KT 15 GeV , 

T D , = 7.37 x 10~ 17 GeV , 

r D / = 9.98 x 10" 17 GeV . (35) 

We found that the semileptonic decay into jl — l + /2 is strongly suppressed. The decay 
widths to it are about one order of magnitude smaller than ones into je = 3 + /2. This 
predicted suppression needs to be checked in the future experiment. 

In Table || we present our results for the branching ratios of B semileptonic decays to 
lowest lying positive parity charmed mesons. We have taken tb = 1.62 ps. The ratio of 
the two semileptonic rates for B decays into D\ and _D 2 mesons is 



R = W^D^ = • { ) 



We are now at the position to compare the theoretical predictions with the available exper- 
imental data. Recently the ALEPH and CLEO collaborations reported the measurement 
of the branching ratios for semileptonic B decay to excited D mesons. For B —>■ D\lv the 
average value of their results is 

B{B -> D 1 lis) = (6.5 ± 2.0) x 10~ 3 . (37) 

The B — ► D\iv branching ratio has not yet been determined. CLEO set the limit 
B(B -> D*£u) < 0.8% f! 7 !, while ALEPH found B(B -> D* lv) < 0.2% p . 



Our prediction for R is larger than available experimental data though we agree with 
the upper CLEO limit for D\. This leads one to consider 0(l/m c ^) corrections. One may 
consider the mixing of the D\ with the jg = l + /2 meson, but it would worsen our prediction 
since, due to the very small decay amplitudes into the je = l + /2 states, it would lessen our 



10 



prediction for B(B — > D\lv). We are thus lead to consider direct 0(l/m Ci &) corrections 
in the decay amplitudes. A detailed discussion of 0(l/m Ci &) corrections can be found in 



Refs. [H, which indeed enhance the rate to B — > D\lv and lead to the expection that its 



branching ratio is greater than that for B — > D\iv. 

V. CONCLUSION 

In this work we have presented the investigation for semileptonic B decays into excited 
charmed mesons. Within the framework of HQET we have evaluated the universal Isgur- 
Wise functions r(y) and ((y) by using QCD sum rules in the leading order of a s . The r(y) 
and ((y) functions can be well fitted by linear approximations r(y)(((y)) = r(l)(C(l))(l — 
Pr((){y ~ !))• The values of r and ( at zero recoil have been given. From a comparison 
with the results of other approaches we found a broad agreement with some of the quark 
model results for r(l), whereas a smaller value for C(l)- 

We have computed, for the decays B — > (Dx.Dl^tv and B — > (D' , D[)£p, the differ- 
ential decay widths and the branching ratios. Our predictions for D\ are smaller than 
experimental data, even take into account the wide spreading and large uncertainty of 
experimental results. Our predictions for D\ are below the experimental upper bounds of 
CLEO's results. We also predict the relation B(B — > D* 2 ii>)jB{B — > Dx£u), which might 
be larger than available experimental limit. The discrepancy may be attributed to a large 
l/m c correction which enhances the transition to D x . We predict tiny branching ratios for 
D'q and D[, but there is no direct experimental check yet. 



After finishing this paper we have learned of a paper by Colangelo et al. [28], in which 
the semileptonic decays B — > (D' , D[) are studied by using a similar approach. In partic- 
ular, the radiative corrections are included in the calculations of Isgur-Wise function £ in 
their work. 
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APPENDIX A: 

We list here the lorentz structures used in the paper. 
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4= \{y 2 - l)g» v + Wv v + (1 - 2y)*/V l/ - (3yvV + (y - 2)v'»v v 
V6 L 

~\ (f™ p {vp ~ r/py) + e"*>« - v' a y)) v a v' p , 

i(l + y)<f*{vP - v'») - i(l + - t/'y) - ±(1 + 

+ ^(1 - j/)*/ W" + ^(1 - ^uVV" - -(1 + y^'W 
2 2 3 

- -(1 + y)t/ + (u'V^ + uVy , 
3 

^ _ v '» ? 

(y- 1)0"" + u'V - u'V , 
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Figure Captions 

Fig. 1. Feynman diagrams contributing to the sum rules for the Isgur- 
Wise form factors. 

Fig. 2. Prediction for the Isgur-Wise form factors. The upper band 
corresponds to the r, the lower one to the £. 
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TABLES 



TABLE I. Parameters of the form factors r, £ 



r(l) 



Pi 



C(i) 



Pi 



Ref. 



0.74 
0.97 
0.71 
0.76 
0.49 
0.54 
1.14 
0.89 
0.94 



0.9 
2.3 
1.5 
0.97 
0.9 
2.8 
1.9 
1.45 
1.50 



0.26 
0.2 
0.8 
0.42 
0.5 
0.62 
0.82 
0.12 
0.45 



0.5 
1.1 
1.0 
0.97 
0.4 
2.8 
1.4 
0.73 
0.83 



This work 






m 



m 



m 

iii, pi 



TABLE II. Branching ratios for semileptonic B decays to P-wave charmed mesons 

Decay mode Branching ratio Exp. result 

B^D x lv 0.34 0.56 ± 0.13 ± 0.08 ± 0.04 @ 

B -» D\lv 0.52 < 0.8(90% C.L.) Jh 

B -» D' £D 0.019 
B -» D^p 0.025 
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FIGURES 




Fig. 2 
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